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Abstract 

It is shown that the chiral symmetry is spontaneously broken by a constant magnetic 
field in QED. The dynamical mass of fermions (energy gap in the fermion spectrum) 

, r 1/2] 

is rridyn — Cve-Bexp — (it /a) , where B is the magnetic field, the constant C is of 
order one and a = e 2 /Air is the renormalized coupling constant. Possible applications 
of this effect are discussed. 
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The dynamics of fermions in an external constant magnetic field in QED was considered 
by Schwinger long ago Wi. In those classical works, while the interaction with the exter- 
nal field was considered in all orders in the coupling constant, the quantum dynamics was 
treated perturbatively. There is no dynamical chiral symmetry breaking in QED in this 
approximation ||. In this paper we reconsider this problem, treating the QED dynamics 
nonperturbatively. We will show that a constant magnetic field B leads to dynamical chiral 
symmetry breaking in massless QED. The dynamical mass of fermions (energy gap in the 
fermion spectrum) is: 

m dyn ~ CVeB exp 

where the constant C is of order one and a = e 2 /4ir is the renormalized coupling constant 
relating to the scale /x = m<^„. 

The essence of this effect is the dimensional reduction D —>■ D — 2 (3 + 1 — > 1 + 1 
in this case) in the dynamics of fermion pairing in a magnetic field, pointed out recently 
in Refs.[]3], |Q. Actually, we will show that, in Euclidean space, the equation describing 
the Nambu-Goldstone (NG) bosons in QED in a magnetic field has the form of a two- 
dimensional Schrodinger equation: 

(-A + m 2 dyn + V( r ))^(v) = 0. (2) 

Here \l/(r) is expressed through the Bethe-Salpeter (BS) function of NG bosons, A = 
d 2 /dx 2 + d 2 jdx\ (the magnetic field is in the +X3 direction; X4 = it), and the potential 
V(r) is 

yw =^ exp (CH-^)- r2=xl+xl (3) 

where Ei(x) = — J^,dtexp(—t)/t is the integral exponential function and I = | e-S | 1//2 is 
the magnetic length. 

We emphasize that we work in the conventional, weak coupling, phase of QED. That 
is, the bare coupling a^, relating to the scale \i = A, where A is an ultraviolet cutoff, is 
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assumed to be small, <C 1. Then, because of infrared freedom in QED, interactions 
in the theory are weak at all scales and, as a result, the treatment of the nonperturbative 
dynamics is reliable. 

We will consider possible applications of this effect at the end of the paper. 

The Lagrangian density of massless QED in a magnetic field is: 

C = --F^F, V + I \$, {yfD^} , (4) 
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where the covariant derivative is 



D^dr-ieiAf + AJ, A? = ( 0, -|ar 2) ^ 0) . (5) 



Besides the Dirac index (n), the fermion field carries an additional, flavor, index a = 
1, 2, . . . , N. Then, at iV > 2, the Lagrangian density (Q) is invariant under the chiral 
SUl(N) x SUr(N) x Ul+r(1) transformations ( we will not discuss the dynamics related 
to the anomalous, singlet, current in this paper). Since we consider the weak coupling 
phase of QED, there is no spontaneous chiral symmetry breaking at B = ||. We will show 
that the magnetic field changes the situation dramatically: at B ^ the chiral symmetry 
SUl(N) x SUr(N) breaks down to the SUv(N) = SUr + l(N) and, as a result, there appear 
iV 2 — 1 gapless NG bosons composed of fermions and antifermions. As we will see, the BS 
equation for the NG bosons defines the dynamical mass (energy gap) for fermions. 

The homogeneous BS equation for iV 2 — 1 NG bound states takes the form (for a review 
see Ref.0): 

Xab{x,v;P) = -i J d 4 x 1 d 4 y 1 d 4 X2d 4 y2G A A 1 (x,x 1 )K AlBi;A2 B 2 (xiy 1 ,x 2 y2) 

■ XA 2 B 2 ( x 2,y2;P)G Bl B(yi,y), (6) 

where the BS wave function Xab{ x > V'^ P) = (O|^''0a(^)'0b(2/)|-P; /?),/? = 1, 2, ... , iV 2 — 1, and 
the fermion propagator Gab(x,v) = (®\TipA(x)4>B(y)\0)', the indices A = (na) and B = (mb) 
include both Dirac (n, m) and flavor (a, b) indices. Note that though the external field A e * 1 (g) 



breaks the conventional translation invariance, the total momentum P is a good, conserved, 
quantum number for neutral channels M, in particular, for these NG bosons. Since, as will 
be shown below, the NG bosons are formed in the infrared region, where the QED coupling 
is weak, one can use the BS kernel in leading order in a [0: 

K Al B 1 ;A 2 B 2 {xm,x 2 y 2 ) = -47rza<J ai0fl <y 6akl 7^ in3 7^ laroi P AB/ (jte - x 2 )S(x 1 - x 2 )6(y 1 - y 2 ) 

+ ^ia5 aibl 5 b2a2 ^ nmi Y m2 n 2 v A x i ~ X2)S(xi - yi)S(x 2 - y 2 ), (7) 
where the photon propagator 

/ h h \ 1 



(2tt) 4 J V k 2 J k 2 

(A is a gauge parameter). The first term on the right-hand side of Eq.(|7]) corresponds to 
the ladder approximation. The second (annihilation) term does not contribute to the BS 
equation for massless NG bosons (this follows from the fact that, due to the Ward identities 
for axial currents, the BS equation for the NG bosons can be reduced to the Schwinger-Dyson 
equation for the fermion propagator where there is no contribution of the annihilation term 
0). Therefore we shall omit this term in the following. Then the BS equation takes the 
form: 

Xab{ x ,Vi P ) = -^a J d A x 1 d 4 y 1 S AAl (x,x 1 )5 aia2 ^ im XA 2 B 2 { x i,yii P ) 

S Bl B{yi, y)V,j, u (yi - xi), (9) 

where, since we are working to lowest order in a, the full fermion propagator G AB {x,y) is 
replaced by the propagator S of a free fermion (with the mass m = m^ n ) in a magnetic field 

S AB (x,y) = 5 ab exp {^-{x - yYA e ^(x + y)\ S nm (x-y), (10) 
where the Fourier transform of 5* is 



S(k) = J dsexp 



is (k^ k\ — — m d yn 



o 

(fcV -k 3 ^ + m dyn )(l + 1 1 1 2 tan(eBs)) - k ±7± (l + tan 2 (e5s))j (11) 



(here kj_ = (ki,k 2 ), 7j_ = (71,72))- Using the new variables, the center mass coordinate, 
R = (x + y)/2, and the relative coordinate, r = x — y, equation (||) can be rewritten as 

Xnm{R, r; P) = -lira j d 4 R 1 d 4 r 1 S nni (r-R 1 + r i^ in2 x n2 m 2 (Ri, n; P)Y m2mi 

■ S^r-^-R + R^V^-^expl-ieir + nrAfiR-R,)}. (12) 



Here the function Xnm(R, t;P) is defined from the equation 

x1 B (^y;P) = (0\TMx)My)\P;P) = X^ b eMier^Af(R))x nm (R,r;P) (13) 

where X 13 are N 2 — 1 flavor matrices (tr(A /3 A 7 ) = 25g 7 ; /3/y — 1,2, ... , N 2 — 1). The important 
fact is that the effect of translation symmetry breaking by the external field is factorized 
in the phase factor in Eq. (|13|) and equation (|12|) admits a translation invariant solution, 
X nm (R,r; P) = exp(—iPR)x nm (r;P). Henceforth we will consider the case with P — > 0. 
Then, transforming this equation into momentum space, we get: 

/ x f d 2 q ± d 2 R ± d 2 k±d 2 k\\ ~ ( KPTt . 

Xnm{p) = -^a j j— - ,i exp(-iq ± R)S nni [p\\,p± + eA (R) + — 



ln 1 n 2 Xn 2 m 2 (^T^rtn^mim [P\\, P-L + eA eXt (R) — 

■ V flu (k ll -p h k ± -p ± ~2eA ext (R)), (14) 

where p\\ = (p°,p 3 ), pj_ = (p 1 ,^ 2 ). 

The crucial point for the further analysis will be the assumption that m^yn <C y|eS| 
and that the region mostly responsible for generating the mass is the infrared region with 



k^ra dyn <C y I eS | . As we shall see below, this assumption is self-consistent (see Eq.(P). 
The assumption allows us to replace the propagator S nm in Eq. ( fH] ) by the pole contribution 
of the lowest Landau level (LLL). In order to show this, we recall that the energy spectrum 
of fermions with m = md yn i n a magnetic field is || 



E n (p 3 ) = ± x lm 2 dyn + 2\eB\n + p 2 , n = 0, 1, 2, . . . . (15) 
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(the Landau levels). The propagator S(p) can be decomposed over the Landau level poles 

D n (eB,p) 



S(p) = i exp 



p1 x 00 



\eBl 



S ( lT pl-pl-ml yn -2\eB\n 



(16) 



with 



D n (eB,p) = (p° 7 ° - pV + m dyn ) 



(l-i 7 V)^n 2 



PL 

|eB| 



- (l + z 7 V)^n-i 2 



_p1 

|eB| 



PL 

\eB\ 



where L n (x) are the generalized Laguerre polynomials (L n = L° and L" 1 (a;) = by defini- 



tion). Eq.flTBD implies that at m^ n <C y |e£>|, the LLL with n = dominates in the infrared 
region with p<m dyn . 

The LLL pole dominance essentialy simplifies the analysis. Now 



S(p) ~ zexp(-/ 2 pi) V m ? n (l-^V), 



(17) 



where py = p° 7 ° — p 3 7 3 and p| = (p ) 2 — (p 3 ) 2 , and equation fll4| ) transforms into the following 



one: 



P(PbP±) 



2al 2 



e -' 2p i / d z A ± d'k ± d'k\\e- l ' A Hi -i-rVh 



2l ^2, -l 2 A 2 , 



fey + 
W 1 



dyn 



P (fe|i,k ± ) fjj + m fV (l - i 7 V)^(fe|| - pii, k ± - A ± ), 



fe|| m d«/n 



where 



P(P||>P-L) = (P|| - m dyn)x{p){P\\ - W d2/ n). 

Eq.(18|) implies that p(py,p±) = exp(— / 2 p 2 L )<y9(p||), where <p(p\\) satisfies the equation: 



(18) 



(19) 



V(P||) 



ira 



(2tt)< 



7 | 1,1 dyn 



fe|| + 



7"(1 - ^VjPjl^fey -py 



(20) 
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^Uh-P\\) = I ^±exp(-^I)2V(A;|| - P ||,k ± ). (21) 

Thus the BS equation has been reduced to a two-dimensional integral equation. Of course, 
this fact reflects the two-dimensional character of the dynamics of the LLL in the infrared 
region, that can be explicitly read from Eq. (p~7|) . 

We emphasize that the dimensional reduction in a magnetic field does not affect the 
dynamics of the center of mass of neutral bound states (in particular, these (N 2 — 1) NG 
bosons). Indeed, the reduction 3 + 1— > 1 + 1 in the fermion propagator, in the infrared region, 
reflects the fact that the motion of charged particles is restricted in directions perpendicular 
to the magnetic field. Since there is no such a restriction for the motion of the center of 
mass of neutral particles, their propagator must have a (3 + l)-dimensional form [|rTJ. This 



fact was explicitly shown for neutral bound states in the Nambu-Jona-Lasinio model in a 
magnetic field, in 1/N C expansion B and for neutral excitations in nonrelativistic systems 



12] . Since, besides that, the propagator of massive fermions is nonsingular at small momenta, 
we conclude that the infrared dynamics of the NG modes is soft in the present model. This 
in particular implies that the phenomenon of spontaneous chiral symmetry breaking in this 



model does not contradict to the Mermin-Wagner-Coleman theorem [13| forbidding the 
spontaneous breakdown of continuous symmetries at D — 1 + 1. 

Henceforth we will use Euclidean space with k^ = —ik°, where the total momentum P of 
NG bosons equals zero. In order to define the matrix structure of the wave function (p(p\\) 
of the NG bosons, note that, in a magnetic field, there is the symmetry SO (2) x SO(2) x V, 
where the SO(2) x 5*0(2) is connected with rotations in X1-X2 and X3-X4 planes and V is the 
inversion transformation x% —>■ —X3 (under which a fermion field transforms as ip — ► «7s73'0)- 
This symmetry implies that the function <^G°||) takes the form: 

<KP[|) = 75(^4 + «7i72^ + P\\C + ijij2P\\D) (22) 

where p\\ = p 3 7 3 + p 4 7 4 and A, B, C and D are functions of p"j (7^ are antihermitian in 
Euclidean space). 
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We begin the analysis of the equation fl2"0|) by choosing the Feynman gauge (the general 
covariant gauge will be considered below). Then, 

M ., . °fdxexp(— l 2 x/2) 

v Uh-p\\) = J {h - n? + x 

and, substituting expansion ( f22|) into equation (p0|), we find that B = —A, C = D = 0, i.e., 
<f(p\\) = ^75(1 — ilijz), and the function A satisfies the equation 

a f d 2 kA(k) J dxexp(-xl 2 /2) 

(henceforth we will omit the symbol || from p and k). Introducing the function 



(2,r)2 F + ' 

we get (from Eq.(^)) the two-dimensional Schrodinger equation (0) with the potential 





- (24) 

(Ko is the Bessel function). Since — m 2 dyn plays the role of energy E in this equation and 
V(r) is a negative, i.e. attractive, potential, the problem is reduced to finding the spectrum 
of bound states (with E = —m dyn < 0) of the two-dimensional Schrodinger equation with 
the potential (p^) . For this purpose, we can use some results proved in the literature. 
First, the energy of the lowest level E(a) for the two-dimensional Schrodinger equation is a 
nonanalytic function of the coupling constant a at a = |Q. Second, if the potential V(r) 
were short-range, then m dyn (a) = —E(a) would take the form m 2 dyn ~ exp[— l/(aa)] where 
a is a positive constant |14|. However, the potential is long-range in our case. Indeed, using 
the asymptotic relations for Ei{x) f5|, we get: 

V(r) ~ -, r -> oo; 

7i r l 



V(r) + r^O (25) 



where 7 ~ 0.577 is the Euler constant. To find m 2 , (a), we shall use the integral equation 
(53) at p = 0. Then, as a — > 0, the dominant contribution in the integral on the right hand 



side of Eq.(^) is formed in the infrared region with soft k 2 <m 2 dyn : 



A(0) 



a 

2t7< 



A(0) 



d 2 k 



k 2 + m 



dyn 



dxexp(-y/2) J*,,^ 
l 2 k 2 + y ~ 4vr 1 ' 



In 



m dynl 



i.e. 



m d . 



yn 



CJ\eB \ exp 



(26) 



(27) 



where the constant C = 0(a°) remains undefined in this approximation. Note that this result 
agrees with the analysis of Ref.|fT5| where the analytic properties of E(a) were considered for 
the Schrodinger equations with potentials having the asymptotics V(r) — > 1/r 2 as r —>■ 00. 

Let us now turn to considering the general covariant gauge (|5p. As is known, the ladder 
approximation is not gauge invariant. However, let us show that, because the present effect 
is due to the infrared dynamics in QED, where the coupling constant is small, the leading 
term in h\{md yn )i ^ n ( m dyn) — — (tt/o) 1 / 2 , is the same in all covariant gauges. 

Acting in the same way as before, we find that in the general covariant gauge the wave 
function fl2~2|) takes the form (p(p) = 75(1 — 27172) (A + pC) where the functions A and C 
satisfy the equations: 

d 2 kA{k) °[dx(l - Ax/ 2 /4) exp(-x/ 2 /2) 



A(p) 



a 
2^ 2 



k 2 + 



nr 



dyn 



C(P) 



aX r d 2 kC{k) 



Ati 2 J k 2 + m 2 dyn 



2k 2 - (kp) 



(k — p) 2 + x 

/c 2 (kp)l °fdxexp(— xl 2 /2) 



[(k-p) 2 + x] 



(2f 



(29) 



One can see that the dominant contribution on the right-hand side of Eq.(p8|) (proportional 
to [lnm^ n / 2 ] 2 and formed at small k 2 ) is independent of the gauge parameter A. Therefore 
the leading singularity in ln(m^„), ln(md yn ) = —(n/®) 1 / 2 , is indeed gauge invariant. 

This concludes the derivation of Eqs. ([!]) and @ describing spontaneous chiral symmetry 
breaking by a magnetic field in QED. 
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In conclusion, let us discuss possible applications of this effect. One potential application 



is the interpretation of the results of the GSI heavy-ion scattering experiments [16 in which 



narrow peaks are seen in the energy spectra of emitted e + e pairs. One proposed explanation 



1 17 1 is that a very strong electromagnetic field, created by the heavy ions, induces a phase 
transition in QED to a phase with spontaneous chiral symmetry breaking. The observed 
peaks are due to the decay of positronium-like states in this phase. The catalysis of chiral 
symmetry breaking by a magnetic field in QED, studied in this paper, can serve as a toy 
example of such a phenomenon. In order to get a more realistic model, it would be interesting 



to extend this analysis to non-constant background fields [18|. 

Another, potentialy interesting, application can be connected with the possibility of the 
existence of very strong magnetic fields (B ~ 10 24 Gauss) during the electroweak phase tran- 



sition in the early Universe As the present results suggest, such fields might essentially 
change the character of the electroweak phase transition. 

Yet another application of the effect can be connected with the role of iso- and chro- 
momagnetic backgrounds as models for the QCD vacuum (the Copenhagen vacuum pOfl ). 
Also, as has been suggested recently [ST], isomagnetic fields in the vacuum of electroweak 



left -right models can induce the parity breakdown. Our work suggests that such field con- 
figurations may play the important role in triggering chiral symmetry breaking in QCD and 
those left-right models. 

The work of I.A.Sh. was supported in part by the International Soros Science Education 
Program (ISSEP) through grant No.PSU052143. 

Note added in proof: 

We have just finished a complete (both analytical and numerical) analysis of integral 
equation fl23|). The result for md yn agrees quite well with estimate ( |27| ) and is md yn = 



C J | eB | exp — 7r/2y / 7r/2« , where C is of O(a ). Notice that the ratio of the powers of this 



exponent and that in Eq.fl2~T|) is 7r/2\/2 ~ 1.1. 
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